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Abstract 

Magnetic properties under the external field are investigated in low-carrier two-band systems, 
which may explain the nontrivial phase boundary found in temperature vs. magnetic field dia- 
gram discovered in some materials, such as filled-skutterudite compound CeOs4Sbi2. Analysis is 
made both for the periodic Anderson model with the small-dispersive / band and the simplified 
two parabolic band model in the vicinity of the Fermi level. The magnetic susceptibilities are 
calculated by using the random phase approximation. It is shown that the maximum value of the 
magnetic susceptibility perpendicular to the external field is enhanced and yields the anomalous 
phase boundary. By applying the magnetic field, the phase boundary shifts to higher temperature 
region in the insulating state with a small band gap. On the other hand, the similar phase bound- 
ary also appears in the semi-metallic states, in which the structure of the density of states in the 
vicinity of the Fermi level and the finite temperature effect are essential. 

PACS numbers: 71.10.-w, 71.27.+a 
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I. INTRODUCTION 



The low-carrier concentration systems have attracted much attention because of the ap- 
pearance of unusual magnetic and transport properties. Ce monopnictides, such as CeP and 
CeSb, are typical compounds, which show the very complicated phase diagram under the 
external field.-^ Ce monopnictides are the semi-metals and the carrier concentrations are 
about 0.01 per Ce ion, so that the Fermi surfaces of electrons and holes are very small. 

Recently, the Pr-based filled skutterudite compounds have also been investigated actively, 
because of the discovery of fascinating superconductivities and magnetic orders.- 

The anomalous metal-insulator transition is observed in the Ce-based filled skutterudite 
CeOs4Sbi2, where the charge density wave (CDW) or the spin density wave (SDW) transition 
occurs at 0.8 K with a gap opening of quasiparticle bands at low temperature region.- The 
phase boundary shifts to the higher temperature region with an increase of the external 
magnetic field,-£ which is contrasted to the usual phase boundaries, which shift to the lower 
temperature region by the field. CeOs4Sbi2 may also be considered as a low-carrier system 
with the semi-metallic band structure. 

The resistivity of CeOs4Sb 12 increases as p ~ J 1 - - 56 even above the transition temperature 
when no magnetic field is applied^, so that it seems that this material might be considered 
as an insulator. However, from the well-known relationship between the energy gaps and 
the lattice constants^, CeOs4Sbi2 should belong to the semi-metallic group. Therefore, we 
consider that the basic band structure is a semi-metal, and the anomalous increase of the 
resistivity at low temperature may be due to disorder and/or magnetic impurity. Note that 
disorder can affect the transport properties of small carrier system drastically. This problem 
is to be studied separately in the future. 

Previously, we have investigated this anomalous metal-insulator transition by construct- 
ing the effective low energy model for CeOs4Sbi2- We employed the semi-metallic band 
structure motivated by the LDA band structure calculation^ and discussed magnetic prop- 
erties. When the enhancement of the perpendicular susceptibility appears with an increase 
of the external field, the phase boundary between the metallic state and the antiferromag- 
netic insulating (AFI) state shifts to the high temperature region. Thus, the phase diagram 
observed in CeOs4Sbi2 can be qualitatively reproduced. 

Although magnetic properties under the external field have been investigated for the 
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Kondo insulator,-*^ we discuss not only the insulating systems but the semi-metallic sys- 
tems in the present study, and investigate the generalized condition of the appearance of 
the anomalous phase boundary in detail. A part of Sec. [IT] of the present results has already 
been reported in the letter.— 

This paper is organized as follows. We construct the effective model with tight-binding 
band of CeOs4Sbi2 and discuss magnetic properties under the external field in the next 
section. In Sec. Hilt We introduce the two parabolic band model in order to consider the 
magnetic response in detail and discuss the general conditions to realize the anomalous phase 
transitions. Summary and discussions are given in Sec. [TV] 

II. EFFECTIVE MODEL FOR CEOS4SB12 

In this section, we investigate magnetic properties of a filled-skutterudite compound 
CeOs4Sbi2. It is because this material may have the semi- metallic band structure and 
becomes a typical low-carrier system with the anomalous phase boundary. 

The temperature-field (T — B) phase boundary of the filled skutterudite CeOs 4 Sb 12 shifts 
to the higher temperature region with an increase of the external magnetic field.-^ Although 
such an anomalous phase boundary has been observed in quadrupole ordered states^, the 
magnetic susceptibility measurement^ suggests the crystalline field ground state is T 7 , which 
has no quadratic moment. Note that although the Sbi2 cluster has a Th symmetry, J = 5/2 
f 1 state of the Ce atom in the Sb 12 cluster feels only the Oh part of the crystal field. Some 
experiments suggest that SDW transition is realized in this material at Tn = 0.8 K, below 
which the material changes from a metal to an insulator.- 

For the filled-skutterudite compounds, rare earth ions placed in the center of the Sbi2 
clusters form body-centered cubic (BCC) lattice structure, where the / orbitals on the rare- 
earth ions strongly hybridize the p orbitals on X12 clusters. These p orbitals form wide 
conduction bands. We employ the periodic Anderson model consisted of the dispersive / 
band and the topmost of the conduction bands. The Hamiltonian is as follows: 
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FIG. 1: Non- interacting band structure at low energy region in the absence of the external field, 
where a c = 1.0, at = 0.03, a 2 = 0.07, and Ef = 0.3, and V = 0.15. Fermi level lies at uj = 0. The 
inset shows the low-energy band structure. 



where c ka (f^ a ) and Cko- (fka) are creation and annihilation operators of a conduction (/) 
electron with the pseudo-spin a (expressing T 7 doublet) and momentum k. e k (e k ) represents 
the band dispersion of conduction (/) electrons. For the filled-skutterudite compounds, we 
take the BCC tight-binding bands with the nearest and the next-nearest hoppings, e k = a c tk 
and e k = aft^+Ef, where t k = cos(k x /2) cos(k y /2) cos(k z /2)+a 2 (cos(k x )+cos(k y )+cos(k z )). 
Vk is the hybridization between / and conduction electrons. Although the hybridization 
generally has the momentum dependence, we neglect this k dependence (Vk = V), for 
simplicity. n ko . represents a number operator (n kcr = c kCT c kj(7 , 4 ff = fl a fk,a)- U is 

the Coulomb repulsion between / electrons. The last term of the Hamiltonian represents 
the Zeeman effect and B ext is the applied magnetic field along z direction. (3 is the effective 
magnetic moment, which is 5/7 (=gj( J z )=5/6 x 6/7) for Ce. In this section, a c is taken to 
be unity, for simplicity, so that the values of other parameters and quantities will be given 
relative to a c . 

In the non-interacting case, Hamiltonian Eq. ([1]) is easily diagonalized, whose eigenvalues 
are given by 



Taking account of the trend of the lattice constants of Ce-based skutterudite compounds and 
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energy gaps, CeOs4Sbi 2 should belong to the semi-metallic group.- Figure [T] shows the non- 
interacting band dispersion in the absence of external magnetic field at low-energy regime. 
Although we employ the simple model, the obtained band structure captures the essential 
features in comparison with the result of the LDA calculation^ at low energy region, where 
the lower band at T point and the upper band at H point slightly intersect the Fermi level. 
Our obtained band dispersions also show the semi-metallic character. 

In general, the Coulomb repulsion between / electrons is screened in the rare-earth metals. 
Furthermore, the effects of the Hubbard-type interaction may be suppressed by the external 
field. Therefore, in order to investigate the Hamiltonian Eq. (TjQ), we employ the random 
phase approximation (RPA), which is one of the mean-field approximations and is justified 
in the weak coupling regime. Then, magnetic susceptibilities for x and z directions are given 
by 

^=^=fd^' (3) 

, (n] = ^ 2 X? T (q) + X; i (q) + 2^? T (q)x? i (q) 
[ ) 4 1 - ^X? T (q)xSi(q) ' U 

where xla' * s the band susceptibility and is defined as 

— T 

xL' (q) = -jf J2 si ( k + q> iuj n)gl> (k, w n ) , (5) 

1 

g a (k,iu n ) = ^2 , (6) 

iuj n + n - e f ka - ■ — 

iu; n + fx- e^ a 

where e%. a = e£ + aB cxt and e{ CT = e[ + a(3B cxt . 

Note that the band susceptibility x(q) consists of 4 x 4 matrix for each momentum q in 
this treatment in an exact formulation. Since / electrons are strongly localized in comparison 
with conduction electrons, the /-/ component of the band susceptibility is dominant for the 
effect of the external field. Furthermore, since only the / electrons have the strong Coulomb 
repulsion, the magnetic instability mainly results from / electrons. Hereafter, we consider 
only the /-diagonal component, for simplicity. 

Figure [2] shows the momentum dependence of the band susceptibility. The maximum 
peak appears at q = — (1, 0, 0). Let us discuss the instability for the ordered state. In the 
absence of the external field, the band susceptibility Xaa' * s independent of the spin index. 
Considering Eqs. ([3]) and (j3J) within RPA, the maximum peak at H point indicates that 
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FIG. 2: Static band susceptibility for various momenta in the absence of the external field at 
T = 0.004. 

the ordering vector is Q = ^(1,0,0), at which x ± (Q) an d X 2 (Q) are strongly enhanced 
with an increase of the Coulomb interaction. This behavior can be explained by the nesting 
property. Since the Fermi surfaces become small spheres at the top of the lower band (r 
point) and at the bottom of the upper band (H point), the particle-hole excitations can 
easily occur between these momenta, which can be easily understood in Fig. HJ The recent 
neutron scattering experiment for CeOs4Sb 12 shows that the ordering vector in the absence 
of the external field corresponds to — (1, 0, 0).— Therefore, although the LDA results have a 
conduction band with a concave dispersion near T point in addition to the band we assumed, 
it may not affect magnetic properties of CeOs 4 Sb 12 . 

Next, we investigate the magnetic response under the external field. Figure [3] shows 
the largest susceptibilities for the parallel component (x z ) and perpendicular component 
(x^) to the external field. For the perpendicular component (x ± ) , the maximum occurs at 
Q = ^(1,0,0), which does not change as far as £> ext < 0.2. On the other hand, the peak 
position of the parallel component gradually changes with an increase of the external field. 
The largest value of x is enhanced within small B ey± range because the external field may 
change the number of the states close to the Fermi level, which is discussed in Sec. IIHI in 
detail. This result indicates that the instability for the magnetic order is enhanced by the 
external field. The transition temperature Tn of the metal-AFI transition is determined by 
1 — (Tn) = within RPA. Because of the enhancement of x ± i the phase boundary of 
the metal-AFI transition shifts to the higher temperature region with an increase of B ext , 
which is shown in Fig. HI 
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FIG. 3: Largest magnetic susceptibilities as a function of the external field at T = 0.003 where 
{7 = 0. The solid (dashed) line represents perpendicular (parallel) component to the external field. 



FIG. 4: Temperature-external field phase diagram at U = 0.2 (solid) and U = 0.22 (dashed). The 
inset shows the experimental result.— 

Note that except for B ext = 0, the direction of the magnetic moments are perpendicular 
to the external field in the ordered state in the present study. In the very weak field region, 
the shift of the phase boundary is independent of the external field. Namely, the slope 
of the phase boundary as a function of temperature becomes infinite, which results from 
the symmetry of the Hamiltonian with respect to the reversal of the direction of the field. 
However, since magnetic anisotropies exist in the realistic materials due to the coupling to 
the lattice, the magnetic moment may not suddenly be directed to x or y axis, but only be 
slightly canted from z direction. 

Further increasing the external field, the phase boundary is closed. The obtained anoma- 
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lous T—B phase diagram can quantitatively reproduce the experimental result of CeOs4Sbi2, 
which indicates that the effective model in this section may describe qualitative properties 
of this material. 



III. GENERAL CONDITIONS FOR ANOMALOUS PHASE DIAGRAM 

In the previous section, we have obtained anomalous T—B phase diagram by constructing 
the effective model of CeOs4Sbi2, where the enhancement of the perpendicular susceptibility 
with an increase of the external field was essential. In this section, we investigate in 
detail the general condition to realize the anomalous phase diagram in the two-band system 
with low carriers, such as the semi-metallic system shown in Fig. [TJ For this purpose, we 
introduce a simplified model Hamiltonian. Since the external field dependence of the band 
susceptibility is directly related to the AFI-metal phase boundary within RPA, which was 
discussed in Sec. HI1 hereafter we focus on the band susceptibility. 

The energy dispersions of the semi-metallic system close to the Fermi level are the almost 
parabolic structure with respect to momenta k with low carriers (Fig. [T] inset). Therefore, 
we approximate the effective band dispersion close to the Fermi level by the two parabolic 
bands. Then the Hamiltonian for the band dispersions is given as 



where a kaa . (awr) is creation (annihilation) operator of an electron with the pseudo-spin a, 
momentum k, and band index a{= 1,2). We assume 2m = 1 for simplicity, and 7 stands 
for the curvature of the band 2 relative to that of the band 1 and the bottom of the band 2 
is located at k = Q. A stands for the overlap (A > 0) or gap (A < 0) between two bands. 
Figure 0(a) shows the schematic band dispersions of Eq. (jSJ) in the absence of external field. 

Then, the perpendicular and parallel components of the band susceptibility are defined 

as 
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FIG. 5: (a) Schematic band dispersions in the absence of external field (A > 0) for one dimension, 
(b) and (c) stand for the shifts by the external field in the band off-diagonal components x]f an d 
x{? of Eq. (jlip at A > 0. Bold arrows stand for the shift directions due to the external field. 

xo(q) = E^f(q), (io) 

after 

k e k+qcr e ka' 

Z,aat n \ _ - 1 \ - /fck+gt) - /( e k|) / 19 x 
JV k C k+qT C k| 

where /(e) stands for the Fermi distribution function. The schematic behavior for the 
external field is shown in Fig. [5](b) and (c). Note that the amplitude of the band diagonal 
susceptibility (Eq. (fT2l) ) is small in comparison with that of the band off-diagonal component 
under the weak external field. 
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A. Absolute Zero Temperature 



At absolute zero temperature (T = 0), we can analytically calculate Eqs. (!9|)- (|T2|) . The 
semi-metallic state (A > 0) has small hole (band 1) and electron (band 2) Fermi surfaces. 
In the absence of the external field, the volume of the Fermi surface of the band 1 is equal 
to that of the band 2. The radius of each Fermi surface is kp a = kp a = (A > 0), where 
kp a is Fermi wave vector of band a with the pseudo spin a. Therefore, the perfect nesting 
of the Fermi surfaces is realized, so that the band susceptibility always diverges at q = Q 
regardless of any spatial dimensions, where Q becomes the ordering vector. If the orbital 
degeneracy remains, the perfect nesting condition may vanish even in the absence of the 
external field. The effect of the orbital degeneracy will be discussed in later section. 

The radius of each Fermi sphere is gradually modified as the external field is applied. 
When both bands intersect the Fermi level which satisfies the following relations A/2 + 
aB cxt > and A/2 — aB cxt /'-f > 0, both Fermi wave vectors at small B cyd are given by 

kk = V j + °B ext 

7,2 _ /A crBcxt" 
Fcr V 2 7 ' 

If a larger external field is applied, band 1 with pseudospin j and band 2 with pseudospin 
| can not intersect the Fermi level, where the Fermi surfaces of band 1 with pseudospin j 
and of band 2 with pseudospin f are lost. 

If the curvature of the band 2 is equal to that of the band 1 (7 = 1), the volume of Fermi 
sphere of the band 1 with spin a is always equal to that of the band 2 with spin a even in 
the presence of the external field, where the perfect nesting still remains between different 
bands. Therefore, when the infinitesimal Coulomb interaction is introduced, the magnetic 
susceptibility always diverges regardless of the presence of the external field within RPA. 
Since the perfect nesting vanishes for 7 ^ 1 and B ext 7^ 0, the band susceptibility remains 
finite, so that magnetic instabilities are suppressed in comparison with the case of 7 = 1. 

On the other hand, if A < where the band insulator with gap |A| is realized, the 
amplitude of the band susceptibility becomes small. When the gap vanishes by applying the 
external field, the band susceptibility Xo(Q) diverges (7 = 1) and becomes enhanced to a 
large value (7 ^ 1). 

However, our interest is the origin of the anomalous behavior of temperature-field phase 
diagram. In next subsection, we discuss the magnetic response to the external field at finite 
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temperature. 

B. Low Temperature Susceptibility 

Magnetic order cannot occur at finite temperature except for the three-dimensional case. 
However, in order to clarify the origin of the anomalous phase boundary, the magnetic 
susceptibilities for various dimensions are investigated in detail within the mean-field level 
in the very low-temperature region (T = 0.01) in this subsection. 

Since, unfortunately, it is not easy to carry out the analytical calculations of the band 
susceptibility at finite temperature, we numerically evaluate Eqs. (l9l- (fT2l) . Note that since 
we employ the parabolic band in energy dispersions, the bandwidth becomes infinite. In 
order to maintain the finite bandwidth, we restrict the wave- vectors within \k x \ + \k y \ < n, 
\k y \ + \k z \ < it, and \k z \ + \k x \ < n in numerical calculations. However, this restriction is 
irrelevant for magnetic properties at low temperature and weak external field. 

The density of states (DOS) reproduces p{u) ~ la*!^ 2 " 1 at uj ~ for d dimensions at 
A = 5 cxt = 0, which is the well-known result. Because of the restriction of the wave- vectors, 
DOS rapidly decreases at sufficiently large \u\ (> 5) region even in two and three dimensions. 

First, let us compare the amplitude of the perpendicular susceptibility Xo{°l) with that 
of the parallel susceptibility Xo(<l) under the magnetic field for various spatial dimensions 
at 7 = 1. Figure [6] shows the largest values of the band susceptibilities at T = 0.01 in 
the semi-metallic case for various dimensions. The perpendicular components of the band 
susceptibilities become always larger than those of the parallel components in whole range 
of the external field and all dimensions. The largest values of the perpendicular components 
for all dimensions appear at q = Q, which corresponds to the nesting vector, while the 
ordering wave vectors of the parallel components are gradually shifted as the external field 
increases. For the insulating case (A < 0), the perpendicular susceptibilities are always 
larger than the parallel susceptibilities. 

Since these behaviors are similar to the results of Sec. [Til qualitative properties of the 
anomalous phase-diagram obtained in Sec. [TTJmay be captured even in this simplified model. 
The obtained result indicates that the magnetic instability perpendicular to the external field 
occurs in this system. Therefore, we hereafter focus on the behavior of the perpendicular 
susceptibility Xo (q)- 
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FIG. 6: Largest values of perpendicular and parallel components of band susceptibilities as a 
function of external field for (a) one dimension, (b) two dimensions, and (c) three dimensions at 
T = 0.01, A = 1, and 7 = 1. The solid (dotted) lines stand for perpendicular (parallel) components. 
The wave vectors in all perpendicular components are q = Q. 

Although the perpendicular susceptibility ^(q) consists of the band diagonal and the off- 
diagonal parts, the amplitude of the latter is greater than that of the former. Thus, magnetic 
properties are almost determined by the off-diagonal part, which consists of two components, 
x]f and x{? terms. The schematic behavior of the band off-diagonal susceptibilities is shown 
in Fig. E(b) and (c), where there are contributions from the band dispersions shifting away 
from the Fermi level and approaching the Fermi level with an increase of the external field. 
Although the perfect nesting still exists at 7 = 1 even under the external field, the amplitude 
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FIG. 7: Largest values of perpendicular susceptibility Xo an d band off-diagonal components as a 
function of external field at A = 1, T = 0.01, and 7 = 1 at low field region. The solid, dashed, 
and dotted lines stand for Xo (total), %i| (Fig. [3(b)), and Xn (Fig- E( c ))> respectively. All wave 
vectors are q = Q. 

of the band susceptibility becomes finite due to the finite temperature effect even at q = Q. 

In order to understand properties of the perpendicular susceptibilities, we show two off- 
diagonal components of the perpendicular susceptibilities in Fig. for the semi-metallic 
case (A = 1). For one-dimensional system, x}\ (Fig- 151(b)) decreases while x{| (Fig- G3( c )) 
increases as the small external field is applied. For two-dimensional system, x\\ ls almost 
independent of the external field, while xj| (Fig. 0(c)) is strongly suppressed at B ext > A/2. 
On the other hand, x]f increases while x\f decreases in three dimensions as the external field 
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is applied. In addition to these behaviors, the conspicuous changes appear at B ext ~ A/2 for 
all dimensional cases. This result indicates that the number of the states close to the Fermi 
level is directly related to the amplitude of susceptibility. With an increase of the external 
field, the number of the states close to the Fermi level for the cases Fig. E(b) and (c) is 
modified. Thus, the weight of particle-hole excitations at low temperature also drastically 
changed by the external field, which affects the amplitude of x\\ an d that of 

For one-dimensional system, DOS p a {w) with band index a has a maximum peak on each 
band edge (u = ±A/2). Therefore, when the overlap between two bands for xj? (Fig. E](c)) 
disappears for by applying B cxt ~ A/2, the amplitude of x\i becomes maximum. For Xn 
(Fig. 0(b)), the change of the number of the states becomes small in comparison with that 
of Fig. EJc). The total susceptibility Xo is mainly determined by the sum of x\1 an d x\p 
so that a enhancement of Xo f°r the external field appears. 

For larger dimensions, the DOS is related to the structure of the perpendicular sus- 
ceptibility Xo- 111 the two-dimensional system, since DOS is independent of frequency u, 
Xu almost becomes constant. However, the overlap between two bands for Xn disappears 
at B cxt > A/2, so that the amplitude is reduced at B cxt > A/2. Therefore the total 
susceptibility xt ls a ^ so reduced at B cxt > A/2. In the three-dimensional case, DOS for 
each band is p\{uj) = \u — A/2 — aB cxt \^ and p^(u) = \uo + A/2 — aB cxt \^. The ratio of 
change of the number of the states to the external field is approximately represented as 
dpM/d^Uo ~ |A/2 + B cxt \-h for Fig. Eflb) and |A/2 - B cxt \~h for Fig. He) at suffi- 
ciently low temperature, so that the decrease of the number of the states for Xn becomes 
slightly larger than the increase of the number of the states for Xty m the small external 
field region. Thus, the total perpendicular susceptibility Xo 1S reduced at B ext < A/2 while 
it is enhanced with further increase of B cxt (> A/2). 

Although the maximum of the band susceptibilities appears at the nesting vector at finite 
temperature, the amplitude strongly depends on the structure of DOS close to Fermi level. 
In the present study, our main objective is to explore the possibility of the anomalous phase 
transition for three-dimensional realistic systems. Hereafter, we discuss magnetic properties 
more in detail in three dimensions. 

In Fig. [8], we show the perpendicular components of the band susceptibilities for various 
A in the three-dimensional system. The cusps appear around at B ext ~ 5, which arise from 
the band edge due to the introductions of the cutoff wave vectors. If the cutoff wave vectors 
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FIG. 8: Largest values of perpendicular components of band susceptibilities Xo as a function of 
external field at T = 0.01 and 7 = 1 for three-dimensional system. The solid, dashed, and dotted 
lines stand for A = 1 (semi-metal), A = 0, and A = — 1 (insulator), respectively. All wave vectors 
are q = Q. The inset shows the low- field part. 

are absent, the amplitude gradually increases for the three-dimensional system in the present 
model. 

However, our aim is to clarify magnetic properties in the low field region. The difference 
in each perpendicular susceptibility particularly appears at B cxt < |A|/2. For the A < 
case, which corresponds to a insulating state with a gap |A|, the enhancements of the 
susceptibilities for the external field can be realized at B ext > |A|/2. When B ext = 0, there 
is no states close to the Fermi level at low temperature due to the existence of the insulating 
gap. When the gap in DOS vanishes due to the introduction of the external field (B exb ~ 
I A |/2), x]l is strongly enhanced. Since there is no states close to the Fermi level for in the 
whole field region, the amplitude becomes almost constant at low temperature. Therefore, 
as far as the insulating state (A < 0) with a sufficiently small gap, the enhancement of the 
magnetic susceptibility Xo can be realized in the low field region. On the other hand, for 
A > 0, each band intersects the Fermi level and the semi-metallic state is realized, where 
the overlap amplitude between two bands is A at B cxt = 0. As the explanation for A > 
is given in Fig. [3, the perpendicular susceptibility xt f°r the external field is suppressed 
within .Bext < A/2. 

Note here that the suppression of the perpendicular susceptibility in the three-dimensional 
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FIG. 9: Largest values of perpendicular susceptibilities x.o as a function of external field for various 
temperature at 7 = 1 at (a) A = —1 (insulator) and (b) A = 1 (semi metal) for three-dimensional 
system. All wave vectors are q = Q. 

system is weakened in comparison with that of the two-dimensional system in the semi- 
metallic states. It is because that the change of the number of the states close to the Fermi 
level becomes gradual. 

C. Susceptibilities at finite Temperatures and Phase Diagram 

Let us consider the possibility of anomalous phase diagrams in the insulating and semi- 
metallic cases in three dimensions. Figure [9] shows the temperature dependence of the 
perpendicular susceptibilities. For the insulating state, the enhancements of the perpendic- 
ular susceptibilities for the external field appear in all temperature regions at B cxt > |A|/2. 
On the other hand, in the semi-metallic state, the perpendicular susceptibility is sensitive 
to temperature. At low temperature (T = 0.01), the decrease of the perpendicular suscep- 
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FIG. 10: Largest values of perpendicular susceptibilities Xo as a function of external field for 
various curvature 7 at A = 1 and T = 0.1 for three-dimensional system. All wave vectors are 

q = Q 

tibility for the external field is clearly shown at B ext < A/2. However, the suppressions of 
the perpendicular susceptibilities are weakened with an increase of the temperature. The 
amplitude is almost independent of B ext (< A/2) at T ~ 0.1. At T > 0.1 the enhancement 
appears in overall range of the external field. The magnetic instability is determined by 
1 — Uxo — * within RPA as discussed in the previous section. If the magnetic order caused 
by the Coulomb interaction is realized at the higher temperature range which exhibits the 
enhancement of the perpendicular susceptibility for the external field, the anomalous T — B 
phase diagram whose phase boundary monotonically shifts to the higher temperature region 
is obtained. 

The effect of the curvature 7 of parabolic dispersions is shown in Fig. (TTJ For all 7, the 
enhancements appear at the low field region. However, the amplitude of the enhancement 
is the most conspicuous at 7 = 1 where the nesting condition survives under the presence 
of the external field. Therefore, as the curvature of one band becomes close to that of the 
other, it is favorable to realize the enhancement of the susceptibility for the external field. 

Finally, let us consider the structure of the phase boundary by taking account of the struc- 
tures of the perpendicular susceptibilities. Figure [TT] shows the possible phase boundaries 
within RPA. In the insulating state (A < 0), the monotonic enhancement of the perpen- 
dicular susceptibility appears within low field region at the wide temperature region. Then, 



17 




FIG. 11: Possible phase boundaries; (a) A < (insulator) or A > (semimetal) and higher 
temperature region, (b) A > (semimetal) and very low temperature region (T <C A). 

obtained phase boundary always shifts to the high temperature region with an increase of 
the external field (Fig. ITTIfa)). For the semi-metallic case (A > 0) in the higher temperature 
region, the structure of the perpendicular susceptibility is similar to the insulating case, so 
that the phase boundary shows also monotonic shift to the higher temperature region for 
the external field. Thus the phase boundary also becomes the structure of Fig. ITTIfa). 

On the other hand, for the semi-metallic case (A > 0) in very low temperature region, 
the perpendicular susceptibility becomes the concave structure within the low field region. 
Thus, when the magnetic order is realized in this temperature region, the phase boundary 
may become nontrivial structure (Fig. ITTT b)). 

To summarize the parabolic band model at finite temperatures, the enhancement of the 
perpendicular susceptibilities by the external field can appear in the insulating states with 
a small gap. On the other hand, for the semi-metallic states when the overlap between 
two bands is sufficiently small, the enhancement may appear except for very low tempera- 
ture region (T C A), In very low temperature region, the perpendicular susceptibility for 
the external field becomes concave structure, so that the phase-boundary may become the 
nontrivial structure. 

In order to realize anomalous phase diagrams in the realistic semi-metallic/insulating 
materials, following conditions are required: (1) both the overlap and gap between two 
bands are sufficiently small; (2) the Fermi surfaces for electrons and holes in the semi- 
metallic state are very small spheres. (3) the change of the density of states close to the 
Fermi level caused by B ext becomes conspicuous; (4) the curvature of one band dispersion 
is close to that of the other. 
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IV. SUMMARY AND DISCUSSIONS 



We have investigated the conditions for the appearance of the anomalous temperature vs. 
external field phase diagram in low-carrier two-band systems with a narrow gap or a small- 
overlap (semi- metal). Employing the periodic Anderson model and the parabolic dispersion 
model close to the Fermi level, the magnetic susceptibilities are calculated by using RPA 
and magnetic properties are discussed under the external field. 

For the filled-skutterudite compound CeOs4Sbi2 whose energy dispersions may have the 
semi-metallic band structure, we construct the effective model by employing the periodic 
Anderson model with the dispersive / band, which reproduces the main features of the LDA 
result close to the Fermi level. Considering magnetic properties, we show that the enhance- 
ment of the perpendicular susceptibility for the external field is related to the anomalous 
phase diagram observed in experimental results. 

We further introduce the two-band model with the parabolic dispersion in order to clarify 
the origin of the anomalous phase boundary. The structure of the perpendicular suscepti- 
bility for the external field is discussed in detail. When a small insulating gap exists in the 
band structure in the absence of the external field, the magnetic susceptibility perpendicular 
to the external field is enhanced by the small field, which results from the disappearance of 
the insulating gap due to the external field. On the other hand, in the semi-metallic case, 
the realization of anomalous phase diagrams strongly depends on the structure of density of 
states close to the Fermi level and the temperature region. 

We consider the effects of Coulomb interaction only by using RPA. This method is one 
of the mean-field approximations and neglects the effect of quantum fluctuations which 
may be particularly important in the vicinity of the quantum critical points.— Therefore, 
if quantum fluctuations are taken into account by employing more accurate methods, such 
as the fluctuation exchange (FLEX) approximation, the phase boundary may quantitatively 
shift to the low temperature region. 

Let us consider the effect of the orbital degeneracy. Heavy fermion compounds have the 
degeneracy of / orbitals of rare earth and actinide ions. We consider the simplest model 
with a lower band and an upper band, where the effect of degeneracy of each band is 
not considered. If one is single band and the others have degeneracy, the perfect nesting 
condition is absent even at zero field. However, in the present study, we consider low-carrier 
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two-band models with a small overlap and a gap between two bands, where the electron 
and hole surfaces become very small spheres. Therefore, when the Fermi spheres become 
obscure at finite temperature region, the nesting vector does not change even if the orbital 
degeneracy is considered. The structure of DOS close to Fermi level becomes most essential 
to the anomalous phase diagram. When the hole doping is also introduced, the nesting 
condition becomes worse. However, as far as the concentration does not conspicuously 
affect the number of the states near the Fermi level, the doping effect is irrelevant at finite 
temperature region in this study. Thus, we conclude that the orbital degeneracy does not 
change the result presented here. 
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